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ABSTRACT 

The study of young massive clusters can provide key information for the formation 
of globular clusters, as they are often considered analogues. A currently unanswered 
question in this field is how long these massive clusters remain embedded in their natal 
gas, with important implications for the formation of multiple populations that have 
been used to explain phenomena observed in globular clusters. We present an analysis 
of ages and masses of the young massive cluster population of M83. Through visual 
inspection of the clusters, and comparison of their SEDs and position in colour-colour 
space, the clusters are all exposed (no longer embedded) by < 4 Myr, most likely less, 
indicating that current proposed age spreads within older clusters are unlikely. We also 
present several methods of constraining the ages of very young massive clusters. This 
can often be difficult using SED fitting due to a lack of information to disentangle age- 
extinction degeneracies and possible inaccurate assumptions in the models used for 
the fitting. The individual morphology of the Ha around each cluster has a significant 
effect on the measured fluxes, which contributes to inaccuracies in the age estimates 
for clusters younger than 10 Myr using SED fitting. This is due to model uncertainties 
and aperture effects. Our methods to help constrain ages of young clusters include 
using the near-infrared and spectral features, such as Wolf-Rayet stars. 

Key words: 


1 INTRODUCTION 


The majority of stars in our universe do not form alone, but 


rather in groups, or clusters (Lada & Lada 20031. There¬ 
fore, it is necessary for us to gain an understanding of how 
these clusters have formed and evolved. The question ad¬ 
dressed in this paper concerns the early phases of cluster 
evolution, specifically how long the clusters remain embed¬ 
ded. The stars belonging to the cluster will form from a mass 
of cold gas, but how long after this initial stage does the clus¬ 
ter keep leftover gas from the formation process bound to 
it? 

To investigate this question we study the young massive 
cluster population of M83. These young massive clusters, or 
YMCs, are often considered to be the precursors to globular 


clusters (GCs) (e.g. Schweizer & Seitzer (19981 and Krui- 
jssen et al. (|2015|)), once only thought to have formed in the 


early universe. Considering these objects as analogues, we 
can study YMC formation to understand how CCs formed 
and how they evolved to the state we see them today. 

The host galaxy of our clusters, M83, is a nearby quies¬ 
cent spiral at 4.5 Mpe (Thimet al.|2003l, with a rich cluster 
population that has been widely studied, for example, |Bas 


tian et al. (2012), Chandar et al. (20101 and Whitmore et al. 


(20111 are only a few recent papers on this topic. M83 has 


been widely studied partly due to this rich population of 
clusters (7290 clusters and possible associations within the 
catalogue used in this paper from Silva-Villa et al. (2014|), 
which means in depth studies can be carried out and popu¬ 
lation statistics are meaningful. There are also clusters of a 
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wide range of masses and ages. Additionally, M83 is thought 
to be very similar to our own Milky Way galaxy, and so any 
results regarding the cluster population of M83 can be used 
to infer results for the Milky Way, whose complete cluster 
population cannot be determined or studied due to the clus¬ 
ters’ and our position in the disc | Larsen|[2008 1. 

Studies of other young massive clusters both within and 
outside the Milky Way show the importance of feedback 
within star forming regions and clusters and give indication 
of the timescales for gas to be expelled from clusters. For 
example, the Arches cluster towards the centre of our galaxy 
has an age of 2.5 — 4 Myr | Blum et al. pooI| and is observed 
to be free of gas. Other clusters of similar masses, such as 
NGC 3603 (aged ~ 1 — 2 Myr) or Westerlund 1 (aged ~ 
3 — 6 Myr) ([Negueruela et al. 20101 (Gennaro et al. 20111 
I Brandner et al.|2008 1 also have no dense gas bound to them. 
It has been found that even for lower mass systems than 
those already mentioned, the clusters can potentially remove 
gas in < 1 Myr ( Seale et al.||2012 K However this is not the 
case for all lower mass regions, as some seem to require more 
time to clear gas, such as W3-main (e.g. Bik et al. (2014|) 
or NGG 346 in the SMC (e.g. Smith ( |2008| ). It has also 
been found that feedback can prompt further star formation 
around the central s ource, such as in the Ruby Ring region 
(Adamo et al.|2o'l2l or Cluster 23 in ESO 338 (Ostlin et al. 


2007). 


This study focuses on the Ha morphology around the 
YMCs. Ha emission traces the gas associated with the clus¬ 
ters, so through visual inspection of how the Ha is orien¬ 
tated, it is possible to determine whether gas is still con¬ 
tained within the cluster or has been expelled. By comparing 
the Ha morphology with the age of the clusters we can find 
the age by which the clusters have removed any remaining 
gas. 


This technique has previously been used by | Whitmore] 
|et al.| ( |2011[ ), who investigated clusters in Field 1 of M83, 
as shown in Fig. They used Ha morphology as a poten¬ 
tial age dating technique for clusters, rather than specifically 
studying the length of the embedded phase. This study pro¬ 
vided an excellent starting point for Field 1, which we now 
expand into the other 6 fields. 

In § [^we discuss the data used for the study and how 
this was used to identify and classify clusters. Additionally, 
we discuss issues with age dating young clusters and how 
we tried to mitigate this problem. § describes our results 
for the length of the embedded phase for clusters and other 
methods for constraining cluster ages. In §|^we discuss the 
implications of the results for GC formation scenarios, in 
particular the Fast Rotating Massive Star scenario (FRMS) 
while in §[^we give our conclusions. 


2 OBSERVATIONS AND TECHNIQUES 
2.1 Data 

Studying the cluster population of M83 requires sufficiently 
high quality data to resolve individual clusters. In addition, 
there is a further constraint as the detail of the Ha mor¬ 
phology needs to be visible around each cluster. In order to 
satisfy these requirements, HST WFC3 data was used, as 
it provides sufficient resolution at the distance of 4.5 Mpc. 


The images were obtained from the publicly available HST 
Legacy Archive for all seven fields of M83 in the U (F336W), 
B (F438W), V (F547M) (F555W for Fields 1 and 2), Ha 
(F657N), I (F814W) and H (F160W) bands. The archive 
provides fully reduced and drizzled data. The orientation of 
the seven fields is shown in Fig. taken from Adamo et 
al, in prep. No transformation was made in the data to the 
Gousins-Johnson filter system, however we adopt U, B, V, 
Ha, I and H for shorthand, as the WFG3 filters correspond 
approximately to these wavelengths. 


Bastian et al. (2012) and Silva-Villa et al. (2014) used 


these images to identify an extensive catalogue of clusters 
within M83 using a combination of automated and manual 
procedures. Photometry in the U, B, V, Ha and I bands was 
already available and we added H band magnitudes for this 
study. Photometry was performed for each of the existing 
clusters on WFG3 IR channel H band images, and added 
to the catalogue. The coordinates were transposed to the 
H band images using the xyxy transform in the astronomy 
library for IDL, and DAOphot was used to find the mag¬ 
nitudes with a 1.5 pixel aperture and 2.4 pixel background 
apertures with annuli of 0.9 pixels. The pixel size of the IR 
detector is 0.13”, giving aperture sizes of ~ 0.2”, or ~ 4.2 
pc. We chose this small aperture size so that they were phys¬ 
ically the same size as the 5 pixel aperture used for the other 
bands, ensuring that there would be no effect from spatial 
differences. Some of the sources were not bright enough in 
the H band to calculate magnitudes and the H band im¬ 
ages covered a slightly smaller area of sky than the other 
Hlters, so several clusters outside the limit of the H band 
had unreliable photometry. Despite these issues, ~ 5600 of 
the 7290 clusters in the catalogue were updated with H band 
magnitudes. 

After removing all sources with poor or no photometry, 
an aperture correction of 0.62 (calculated from a sample of 
approximately 20 sources) was subtracted to correct for lost 
flux due to a small aperture. Furthermore, Fields 1 and 2 
were dithered differently to the other 5 fields. Our very small 
apertures were highly susceptible to this difference, resulting 
in a fairly significant reduction to the magnitudes in Fields 3 
to 7. We corrected for this by subtracting 0.6 from the mag¬ 
nitudes in those helds. This value was found from comparing 
magnitudes of the same sources in the overlapping regions 
between Fields 2 and 4. The H band photometry is used as 
a method for age dating clusters, as described in § |2.4.4[ 


2.2 Cluster selection 


As part of the selection process [Bastian et al.| ( |2012[ ) and 
Silva-Villa et al. (2014) used concentration indices (GIs) to 


rehne the sample. A GI is a measurement of how centrally 
concentrated emission is for a source. This value helps dis¬ 
tinguish between field stars and clusters; stars should be 
highly centrally concentrated as point sources, whereas clus¬ 
ters have extended emission. Star contaminants can there¬ 
fore be removed by cutting sources below an index selected 
from histogram plots. A histogram plot of the concentration 
index for each source demonstrates two defined sections: the 
hrst is a sharp peak consisting of individual stars followed 
by a lower peak with a smooth tail that illustrates the clus¬ 
ter population. Removing all sources with values in the first 
peak leaves the cluster candidates. These candidates were 
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Figure 1. The orientation of the seven MSS fields in HST WFC3. 
As can be seen, the majority of the galaxy is covered, with a few 
outer regions missed. Image taken from Adamo et al. (in prep). 


then visually inspected to determine if the source was a re¬ 
liable cluster (class 1), potentially a clnster (class 2) or an 
incorrectly identified as a cluster in the previous steps (class 
3). The last of these populations were removed. Error cuts 
were also applied to remove sources with poor photometry 
and aperture corrections were applied to account for flux 
missed by the aperture. We used the most recent updated 
version of this catalogue from Silva-Villa et al.| | 20141, which 
included photometry in the U, B, V, Ha and I bands as well 
as estimates of the age, mass and extinction {Av) for each 
cluster. These properties were obtained by comparing each 
cluster’s spectral energy distribution (SED) with models of 
simple stellar populations as per [Adamo et al.| ( |2010a,b] ) . 

To identify the possible YMCs that could be considered 
GC progenitors, age and mass cuts were applied to the cat¬ 
alogue. No distinction was made between class 1 and 2; all 
were included in our final selection. Our final sample con¬ 
sisted of 91 clusters that were younger than 10 Myr and 
had a mass > 5000 M0. The lower mass cut was applied to 
minimise the effect of stochastic sampling. 


2.3 Cluster categorisation 


After identifying the clusters, we determined whether each 
would be classed as embedded (no shells or areas where gas 
had been cleared), partially embedded (part of the cluster 
had cleared gas) or exposed (clear bubble around the cluster, 
possibly surrounded by a shell of ionised gas). Categorisation 
of the clusters was carried out by visual inspection of the 
Ha morphology associated with each cluster. Images of each 
cluster were created by cropping the original fits files to 
a 12 x 12 arcsecond area with the cluster at the centre. 
After creating images in each band, the B, Ha and I bands 
were stacked to create colour images for each cluster, clearly 
displaying the position of gas in the cluster (traced by the 
Ha emission). Ha continuum subtracted images (provided 
by William Blair, and as discussed in Blair et al. (2014l) 
were used to confirm whether the cluster was embedded or 
not, if the colour image was unclear. 

Examples of each category are shown in Fig.|^ The first 



Figure 2. Examples of clusters in different categories. The num¬ 
ber in the right hand corner gives the cluster’s ID within the cat¬ 
alogue and the concentric circles indicate radial distances from 
the centre of the cluster to 10, 20, 40 and 80 pc respectively, to 
indicate how far the cluster has cleared gas. The top image is 
an embedded cluster, the second cluster is partially embedded 
and the bottom cluster exposed, shown by the clear bubble. Blue 
colour traces B band emission, green traces Ho and red I band. 


cluster is classed as embedded, as it is clearly surrounded by 
gas (represented in green), which is reinforced by a high Ay 
of 2.70, caused by dust associated with the gas. SED fit¬ 
ting gives the cluster a very young age of ~ 4 Myr, and 
therefore it is reasonable to expect that the cluster is still 
embedded, although we will return to issues of age fitting in 
§ |2.4| The second cluster is partially embedded, with a lower 
Av of 0.39 and an age estimate of ~ 7 Myr. These interme¬ 
diate clusters provide a snapshot of the way in which the 
clusters have removed their gas to progress to an exposed 
phase. They can be used to provide an estimate of the age 
at which clusters begin to expel their excess gas and indi¬ 
cate that clusters above this age should be at least partially 
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Figure 3. Illustration of the process used to obtain an indica¬ 
tion of the extent of clearing of the gas from the clusters. The 
two shaded areas show the angles that were measured for the 
minimum and maximum distances cleared in red and blue re¬ 
spectively, labelled with the distance. Angles and pixel lengths 
were measured in GNU Gimp. 


free of gas. For this reason, they could be considered as the 
most important clusters in the sample, however the classi¬ 
fication between partially embedded and either embedded 
or exposed is sometimes difficult to discern, depending on 
how advanced the clearing is for the cluster. A cluster al¬ 
most entirely devoid of gas may appear to be only slightly 
embedded but the gas may not be associated with the clus¬ 
ter at all. The final cluster is exposed with an age of ~ 6 
Myr, as shown clearly by the bubble and shell of ionised 
gas, which has been cleared to a distance of > 80 pc on the 
northern side, but only ~ 30 pc on the southern side. This 
is an example of a ‘blowout’. 

The poor quality of the photometry (offset location on 
the colour-colour diagram possibly caused by contamina¬ 
tion from neighbouring sources) for some clusters produced 
wrong estimates for their physical properties via SED fitting. 
These outlier sources could also be dominated by stochas¬ 
tic IMF sampling, which is known to produce extreme off¬ 
sets with respect to traditional stellar evolutionary tracks 
( Fouesneau fc Langon[ 2010K Therefore we excluded these 
sources. After identification they were confirmed to be un¬ 
reliable sources from visual inspection of the images, where 
all the clusters appeared to have small radii and were red 
in colour. This may have been due to contamination from 
nearby sources. Several further clusters were removed due 
to poor quality images, possibly caused by the continuum 
that is also detected in the Ha band, meaning that the ori¬ 
entation of the Ha could not be accurately determined from 
visual inspection. After these clusters were removed from 
the subsample, the final confirmed YMCs numbered at 35 
exposed, 16 partially embedded and 15 embedded clusters. 

After identifying the state of embedded-ness for each 



Figure 4. The distribution of ages for the young massive clusters 
using SED fitting of U, B, V, Ho and I bands. The orange solid 
line is embedded clusters, the purple dashed line is partially em¬ 
bedded clusters and the teal dot-dashed line are exposed clusters. 


cluster, each of the partial and exposed clusters were studied 
again to find the extent of the gas cleared. Using measure¬ 
ment tools in the GNU Gimp image processing program, the 
distance from the centre of the cluster to the closest point 
on the inner bubble edge was measured in pixels. The dis¬ 
tance to the part of the inner bubble edge furthest from the 
centre was measured in the same way, unless the distance 
exceeded the largest circle on the image, in which case the 
distance was recorded as > 80 pc, for simplicity. Pixel mea¬ 
surements were converted to parsecs to give the distances. 
Finally, as an indication of the extent to which the bubbles 
had been cleared, the angular area of the circular region 
surrounding the cluster that had been cleared to these max¬ 
imum and minimum distances were measured using tools in 
GNU Gimp, as shown in Fig. If there were multiple re¬ 
gions cleared to similar distances, the angles measured for 
each were summed. These angular areas were converted into 
percentages of the area surrounding the cluster by dividing 
by 360'’, as included in Table This information was used 
to compare with masses and determine whether a correla¬ 
tion exists between mass of the cluster and the size of the 
bubbles, as discussed in § |3.3| Additionally, the percent¬ 
age information could potentially indicate the existence of 
a blowout in the gas. We did not measure inner and outer 
shell radii for embedded clusters. 

Two of the partially embedded clusters in Table as 
indicated by *, actually have no maximum distance cleared. 
These two clusters initially look embedded and are sur¬ 
rounded by gas, however, after inspection of the clusters’ 
SEDs and their position on the colour-colour plot (as dis¬ 
cussed later in § |2.4| and § |2.4.1| respectively) it was clear 
that they were not completely obscured by gas, and had 
likely cleared an area in our line of sight. This caused them 
to appear bluer and could be accurately identified as par¬ 
tially embedded. 
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2.4 Refined age estimates of young clusters 

Fig .[^ shows the distributions of SED fitted ages for embed¬ 
ded, partially embedded and exposed clusters. It is impor¬ 
tant to note that the SSP models used to ht the ages include 
nebular emission, which should improve accuracy, and an 
assumed covering fraction of 50% (50% of the ionising pho¬ 
tons are absorbed within the aperture). Additionally these 
ages can be affected by an age/extinction degeneracy and 
model uncertainties due to factors such as covering fraction 
or metallicity, for example. As demonstrated, the embedded 
population peaks at a younger age than the partially embed¬ 
ded and exposed clusters, as we would expect. The partially 
embedded distribution peaks at a more advanced age than 
would be expected, at the same age as the exposed popula¬ 
tion. This could indicate that the process of removing gas is 
very rapid, or just that the earlier stage partially embedded 
clusters appear embedded and have been wrongly identihed. 
Alternatively, this could be an effect of binning the ages into 
the distribution. As this effect is Ha morphology dependent, 
we found that this could be due to an issue with Ha. 

The clusters in this sample have been age dated using 
SED fitting of U, B, V, Ha and I bands. The degeneracies 
associated with this method can give inaccurate estimates 
for the ages of young clusters less than ~ 10 Myr old. While 
using all five of the aforementioned bands for the fitting pro¬ 
vides the best age estimates for the majority of the clusters, 
including Ha presents a problem with some of the partially 
embedded and exposed clusters. 


2 . 4.1 Colour-colour plot 


An important tool in determining the ages of clusters, which 
can also be used to improve age estimates and that displays 
the problems experienced in using SED fitting is the colour- 
colour plot. Fig. shows this plot for the YMC sample of 
M83. The model plotted over the points is the cluster evo¬ 
lutionary model for M83 at solar metallicity that includes 
nebular and continuum emission and uses a Kroupa IMF 


from Zackrisson et al. (20111, which was also used for SED 


fitting of the clusters. The orange circles correspond to em¬ 
bedded clusters, the purple diamonds to partials and green 
stars to exposed clusters. Additionally, we have labelled sev¬ 
eral of the logarithmic ages along the model track. There are 
some key differences between these populations, which are 
clearly visible on the plot, and which can provide further 
information about the early evolution of these clusters. 

The embedded clusters lie primarily in the top right 
hand corner of the distribution. This is expected, as they all 
have young ages so should lie on the lower end of the model, 
but are highly extincted due to their surrounding gas and 
dust and are moved to the top right section. The exposed 
clusters also lie where we would expect, in the bottom left of 
the diagram. This indicates that they are still very young, 
however, unlike the embedded clusters, a lack of gas means 
they have much lower extinction and so are much closer to 
their position on the model. The partially embedded clusters 
lie scattered around the plot, with extinction varying due to 
the orientation of the exposed part of the cluster. If the clus¬ 
ter has cleared our line of sight from gas, then the extinction 
will be lower than a cluster that has cleared in another di¬ 
rection. The position of points on the graph shows a positive 



Figure 5. Colour-colour plot of U-B vs V-I for the final sample 
of YMCs in MSS. The model track by Zackrisson et al, 2011 is 
overplotted to indicate where on the evolutionary track the clus¬ 
ters would lie. Orange circles indicate embedded clusters, purple 
diamonds are partially embedded and green stars are exposed 
clusters. 


result: the classifications of clusters using the colour-colour 
plot agree with those from the morphological analysis. 

The colour-colour plot can also provide information 
about the age that clusters start to remove gas. In the region 
dominated by the exposed clusters, they appear to form a 
line which is approximately parallel to the extinction vec¬ 
tor. This line could be showing the change in the amount of 
gas cleared over time. If the line is followed back along the 
extinction vector to the model track, it indicates that the 
clusters initially began to remove gas at an age of ~ 2 Myr, 
depending on the model. 

This plot also indicates why SED fitting can experience 
issues when fitting ages for embedded and partially embed¬ 
ded clusters. The clusters are brought back along the extinc¬ 
tion vector until they meet the model, giving estimates of 
the age and amount of extinction for the cluster. The shape 
of the model track introduces a potentially problematic de¬ 
generacy for embedded and partially embedded clusters in 
the top right of the distribution. 

The plot indicates that clusters begin to expel gas at 
a few Myr, and it only takes this long to create a shell HII 
region. The embedded clusters have compact HII regions, 
indicating that they should be very young, and therefore 
be fitted ages on the earliest section of the model track, 
however they can be fitted to three different sections of the 
track due to the degeneracy, assigning them lower extinc¬ 
tions and older ages than they actually have. Partially em¬ 
bedded clusters can also be susceptible to this problem. This 
age-extinction degeneracy can only be resolved with accu¬ 
rate knowledge of the UV slope (Calzetti et al.|2014(. 


2 . 4.2 The problem with Ha 

Fig .|^shows an approximation of the Ha excess for the three 
categories of clusters. The Ha excess is the average of the 
magnitudes in V and I minus Ha for the cluster. This value 
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Figure 6. A histogram plot of Ho excess for embedded, partially 
embedded and exposed YMCs in the sample. The orange solid 
line is embedded clusters, the purple dashed line partially em¬ 
bedded clusters and the teal dot-dashed line is exposed clusters. 
They each peak at different levels of excess, with partially em¬ 
bedded clusters split into 2 peaks, indicating differences in the 
Ho morphology. 


Figure 7. New age distribution for the cluster sample, after se¬ 
lecting the correct fit from each cluster’s SED, position on the 
colour-colour plot and Ho morphology. As before, the solid orange 
line is the embedded population, the dashed purple line represent 
the partially embedded sample and the teal dot-dashed line is the 
open clusters. 


provides an indication of the how much the emission from 
the gas has contributed to the Ha magnitude, by subtracting 
the contribution from the stars (approximated by {V -|-/)/2, 
as no other bands nearer Ha were available). Using (U-|-7)/2 
accounts for extinction and offers a close alternative to the 
actual continuum. 

When using SED fitting to approximate ages, young 
clusters would be expected to have a certain amount of 
Ha emission as they may still contain or be surrounded by 
ionised gas. As the figure shows, the exposed clusters peak at 
close to zero excess. This means that even for newly exposed 
clusters with the expelled gas still nearby, no contribution is 
made from the gas to the Ha magnitude. This would suggest 
that these clusters have cleared gas to a radius larger than 
the aperture used for the photometry at a young age. By us¬ 
ing Ha in the fit for the age, the cluster appears older due to 
a lack of flux in this band, despite still being young enough 
to produce radiation that ionises the surrounding gas. In this 
case, the gas is arranged in a bubble surrounding the clus¬ 
ter, and therefore the majority of the Ha emission is from 
a shell outside of the aperture. If the clusters are exposed 
but there is still Ha emission within the aperture from the 
nearby gas, a fit using Ha may still be more accurate than 
without. Bastian et al. (2014a,bl found a similar effect in 
the nebular continuum in the near-infrared when studying 
YMCs. 

The partially embedded clusters have 2 peaks in the 
histogram in Fig.|^ This indicates 2 populations: those that 
have cleared gas in our line of sight and have close to zero 
Ha excess, and those that still have contribution from the 
gas. Each population would require different fits, sometimes 
including Ha and others not. 

The embedded clusters all have Ha excess as the gas 
is still within the cluster, or closely nearby and therefore is 
always included in the measurement of the magnitude. It 


is unlikely that any embedded clusters would require a fit 
without Ha. 

One of the main issues with Ha is the small apertures 
that are used to obtain photometry. They are ideal for cap¬ 
turing the continuum emission, but not the entire emission 
from the region, so the Ha fluxes can be underestimated. 
This loss of flux cannot be accounted for in the aperture 
correction, as the morphology of the gas surrounding the 
clusters that may contribute to the Ha flux varies widely. 
This also means that the excesses displayed in Fig. are 
heavily dependent on the size of the aperture used. Con¬ 
tamination from neighbouring Ha sources can also affect 
the flux in crowded regions. 


2.4-3 Accounting for the problem 

As shown in Table by fitting some of the clusters ex¬ 
cluding Ha contribution (only using U, B, V and I bands), 
their age falls above the cut off for the YMC sample and 
in some cases, the mass falls below the mass limit. Visual 
inspection of the images indicates that, though the optimal 
fit may be without Ha contribution, the new lower masses 
do not seem to visually correspond with the cluster. For ex¬ 
ample, cluster 50574, has a new log(mass) of 2.72, however 
the cluster appears much more massive, and the original fit 
of 4.09 appears more likely. Additionally, cluster 50611 has 
a preferred fit including Ha, however, the image shows no 
Ha emission closely surrounding the cluster, so why includ¬ 
ing Ha is preferred is unclear, though it may be due to the 
continuum. 

It is clear, however, that despite improvement for many 
clusters by tailoring the fit to their morphology, a problem 
still exists with the degeneracy within the colour-colour plot 
and fitting procedures for ages of the clusters. The youngest 
parts of the model track at ~ 2 Myr progress very rapidly. 
This means that, according to the model, clusters are un- 
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likely to be found in this region, and these ages are less 
likely to be assigned to clusters. Additionally the degener¬ 
acy caused by clusters moving through several sections of 
the model track when corrected for extinction mean that 
ages can easily be assigned incorrectly. 

In order to investigate this effect, ages for all the clusters 
in the original catalogue were obtained using SED htting 
both including and excluding Hq flux. The young exposed 
and partially embedded clusters were selected and the 
values for the fit were inspected. The fit that produced the 
minimum was selected as the first optimal fit for the 
cluster. This way of selecting an appropriate ht proved to 
be inaccurate for the majority of the clusters as most of the 
embedded clusters were also assigned non-Ha hts, clearly 
in contradiction with their morphology. These fits also put 
the embedded clusters at a much more advanced age than 
the exposed clusters, which is unphysical. Having found that 
this method was not entirely reliable, a more robust selection 
process was carried out for the fits. A youth indicator that 
is linked directly to the stars is therefore preferred over line 
emission. A potential indicator is the UV band, which will 
be explored further by Calzetti et al. (20141 and the HST 
Legacy Extragalactic UV Survey (LEGUS). 

The position of the cluster in the colour-colour plot, its 
SED and extinction were all used in determining the opti¬ 
mum fit for the cluster. Visual inspection was also used to 
confirm whether the ht seemed reasonable considering the 
Ha morphology and the visually apparent mass of the clus¬ 
ter (higher mass clusters should be noticeably brighter and 
larger in radius, though this is a rough technique used only 
to help conhrm a ht). For example, an exposed cluster may 
be assigned a non-Ha ht due to the yf value, however this 
ht does not account for nearby Ha emission (if the cluster 
is newly exposed and gas is still close), which can be evi¬ 
dent from the strength of the Ha in the cluster’s SED. A 
very large mass may be assigned to the cluster using this 
ht, when the cluster is clearly a lower mass cluster from 
only visual inspection. The SEDs proved the most valuable 
tool in this ht and could even possibly identify the progress 
of a partially embedded cluster. Later stage partials were 
identihed from having very similar SEDs to exposed clus¬ 
ters (much brighter in blue bands) but with slightly higher 
Ha due to nearer gas. Earlier partials looked more akin to 
embedded clusters with high Ha and faint bluer bands, but 
could be distinguished by lower extinction. Several clusters 
were identihed as partially embedded, despite appearing em¬ 
bedded, due to a lack of strong Ha and bright blue bands. 
These seemed to be the clusters that had cleared gas along 
our line of sight but not necessarily in all directions. By em¬ 
ploying all of these methods, a new optimum ht was assigned 
to each cluster, as shown in Table 

Partially embedded clusters are the most difhcult to 
ht. It is often unclear whether a Ha or non-Ha ht is most 
appropriate as they can vary in Ha morphology. This means 
that many of the ages for the partially embedded clusters 
are unreliable and have been affected by the degeneracy in 
the models due to extinction. 

Embedded clusters can also be ahected by degeneracies 
in the htting procedures. Our results show that clusters are 
initially exposed by at least 6 Myr (and most likely much 
younger), and so it is unrealistic for embedded clusters, at an 
earlier evolutionary phase, to be older than exposed clusters, 


as a few are in our initial htting. Due to their high extinc¬ 
tion, they are prone to be assigned older ages due to the 
degeneracy in the model track, meaning there are several 
diherent potential ages for the clusters. In these cases, even 
the Ha hts cannot obtain the younger ages that should likely 
be assigned to these clusters, therefore the ages provided for 
a small sample of the embedded clusters in the sample may 
be too large. 

Using methods like those described previously; consid¬ 
ering SEDs, visual inspection of the image and position on 
the colour-colour plot to determine which ht is appropriate, 
provide a more accurate age estimate for the cluster than 
only SED htting. We identihed 17 exposed clusters that were 
better htted without Ha and were reassigned ages. Most of 
these ages were only marginally different from the cluster’s 
previous age, with only 7 clusters experiencing a fairly sig- 
nihcant change in age. Several of the ages were older, in 
contrast with the apparent age of the cluster from the Ha 
morphology, though the majority were given younger ages. 
In addition to the exposed clusters, three partially embed¬ 
ded clusters were also assigned different ages. 

Fig.[7| shows the new age distribution after selecting 
the correct hts for the clusters. The new distributions peak 
a lot closer together than in the previous plot and there is 
overlap between the diherent populations. The distribution 
for the embedded population of clusters has not changed, 
as they were consistently htted with Ha. Their ages would 
not have improved without Ha, as they would have been 
assigned lower extinction and older ages due to the model 
degeneracy. The partially embedded clusters now have two 
peaks, including a younger sub-population. The latter peak 
consists of those with an optimal ht including Ha, while 
the new younger peak are those without. The exposed clus¬ 
ters now also have several peaks. The oldest dominant peak 
are those still htted with Ha. The two younger peaks align 
with the embedded cluster distribution, indicating a diverse 
range of interaction between the stellar component of the 
cluster with its gas. It could potentially indicate that some 
clusters become exposed at a faster rate than others. The 
overlap between diherent categories of clusters could also be 
another indicator that the transition from embedded to ex¬ 
posed phase is very short, and conhrms that the partially 
embedded clusters represent an intermediate phase of clus¬ 
ter evolution to an exposed state. It also shows that this 
process is continuous and that a majority of the clusters 
have a compact bubble of ionised gas by ~ 3 Myr. 


Infrared age dating of clusters 


Space-based resolution capabilities will soon be limited to 
near-IR data and the traditional extragalactic cluster popu¬ 
lation studies in optical and UV wavelengths will no longer 
be feasible with the retirement of HST. These clusters should 
be visible in the near-IR, however studies that attempted to 
use infrared data, e.g |Grosb0l &: Dottori| ( [2013[ ), using aper¬ 
ture photometry as with optical data, have proven unreli¬ 
able due to spatial resolution effects ( Bastian et al.|[2014a |. 
Despite this, space-based instruments with high spatial res¬ 
olution (e.g. JWST) in the near-IR can still be relevant and 
useful to cluster studies. 

As explained previously, current age dating methods for 
clusters, such as SED fitting, as used for this sample is poor 
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Final sample of young massive clusters in M83 


Cluster ID 

Cat 

SED log(age) 
(Best fit) 

Log(mass) 
(Best fit) 

Av 

Max distance 

(pc) 

Min distance 

(pc) 

% max 

% min 

Best fit 

10141* 

1 

6.60 

4.03 

0 

n/a 

n/a 

n/a 

n/a 

No Ha 

10165 

2 

6.70 

4.00 

0 

> 80 

23.1 

46 

10 

No Ho 

10426 

0 

6.60 

4.00 

2.7 

n/a 

n/a 

n/a 

n/a 

Ho 

10433 

0 

6.60 

4.58 

2.86 

n/a 

n/a 

n/a 

n/a 

Ho 

10438 

1 

6.85 

4.79 

1.52 

28.0 

7.8 

4 

18 

Ho 

10452 

0 

6.70 

4.18 

3.53 

n/a 

n/a 

n/a 

n/a 

Ho 

10594 

2 

6.30 

4.46 

0 

48.9 

14.9 

9 

19 

No Ho 

10597 

2 

6.60 

4.72 

0 

54.9 

10.0 

5 

22 

No Ho 

10645 

2 

6.60 

3.97 

0 

> 80 

10.0 

28 

31 

No Ho 

10732 

2 

6.60 

4.03 

0 

> 80 

20.0 

21 

28 

No Ho 

10915 

1 

6.75 

3.73 

2.78 

10.0 

0 

16 

83 

Ho 

10957 

0 

6.73 

3.57 

2.08 

10.0 

0 

11 

89 

Ho 

11164 

2 

6.60 

4.01 

0 

> 80 

26.5 

16 

28 

No Ho 

11174* 

1 

6.73 

3.95 

0.14 

n/a 

n/a 

n/a 

n/a 

Ho 

11266 

1 

6.88 

3.82 

0.83 

30.9 

0 

17 

21 

Ho 

11297 

1 

6.60 

4.06 

0 

> 80 

0 

14 

13 

No Ho 

11305 

1 

6.70 

3.81 

0.46 

13.8 

0 

8 

92 

Ho 

11312 

0 

6.60 

3.92 

2.7 

n/a 

n/a 

n/a 

n/a 

Ho 

11414 

2 

6.70 

3.87 

0 

> 80 

67.2 

90 

10 

No Ho 

20249 

0 

6.48 

3.98 

3.54 

n/a 

n/a 

n/a 

n/a 

Ho 

20281 

2 

6.60 

4.08 

0 

23.1 

10.0 

8 

17 

No Ho 

20288 

2 

6.78 

4.08 

0 

> 80 

40.0 

45 

14 

No Ho 

20314 

1 

6.73 

3.81 

0 

30.2 

0 

5 

15 

No Ho 

20417 

2 

6.60 

4.26 

0 

> 80 

40.0 

20 

38 

No Ho 

20866 

2 

6.88 

3.75 

0.02 

> 80 

15.6 

33 

35 

Ho 

20868 

2 

6.88 

3.80 

0.18 

> 80 

14.8 

33 

35 

Ho 

20953 

1 

6.73 

3.70 

0.75 

> 80 

53.2 

63 

13 

Ho 

30126 

1 

6.87 

3.97 

0.79 

23.4 

0 

18 

40 

Ho 

30701 

0 

6.60 

3.82 

2.19 

n/a 

n/a 

n/a 

n/a 

Ho 

30989 

0 

6.30 

3.98 

2.29 

n/a 

n/a 

n/a 

n/a 

Ho 

31015 

2 

6.60 

3.93 

0 

> 80 

15.0 

13 

18 

No Ho 

31016 

2 

6.70 

4.00 

0 

> 80 

20.0 

9 

5 

No Ho 

40027 

0 

6.48 

4.19 

3.77 

n/a 

n/a 

n/a 

n/a 

Ho 

40117 

2 

6.87 

4.34 

0.54 

> 80 

15.2 

6 

30 

Ho 

40358 

2 

6.88 

3.89 

1.22 

> 80 

24.0 

32 

7 

Ho 

40543 

2 

7.78 

3.37 

0 

> 80 

46.4 

90 

7 

No Ho 

40820 

2 

6.95 

4.20 

0.38 

> 80 

56.8 

94 

6 

Ho 

40840 

2 

7.30 

3.23 

0 

> 80 

44.9 

78 

7 

No Ho 

40915 

2 

6.78 

3.70 

0 

23.7 

10.0 

27 

11 

Ho 

40926 

2 

6.78 

4.1 

0 

37.5 

7.9 

16 

41 

Ho 

50015 

2 

6.78 

3.88 

0 

40.0 

5.5 

3 

20 

Ho 

50019 

2 

6.79 

3.74 

0 

> 80 

13.3 

45 

7 

Ho 

50035 

2 

6.60 

2.50 

0.78 

> 80 

10.0 

12 

43 

No Ho 

50081 

1 

6.88 

3.77 

0.39 

35.3 

0 

14 

52 

Ho 

50099 

0 

6.60 

3.73 

1.99 

n/a 

n/a 

n/a 

n/a 

Ho 

50140 

2 

6.30 

4.21 

0 

> 80 

43.1 

60 

11 

No Ho 

50160 

1 

6.70 

3.70 

0 

n/a 

n/a 

n/a 

n/a 

Ho 

50233 

2 

6.78 

4.04 

0 

> 80 

7.6 

21 

28 

Ho 

50256 

2 

6.87 

3.81 

0 

> 80 

28.4 

24 

10 

Ho 

50358 

0 

6.70 

3.96 

3.11 

n/a 

n/a 

n/a 

n/a 

Ho 

50377 

1 

6.78 

3.94 

0 

> 80 

0 

11 

22 

Ho 

50574 

2 

6.78 

4.09 

0 

> 80 

45.4 

79 

3 

Ho 

50611 

2 

6.78 

3.70 

0 

> 80 

14.8 

5 

13 

Ho 

50789 

1 

6.60 

4.73 

0 

31.5 

0 

4 

41 

No Ho 

50905 

0 

6.60 

3.70 

0.33 

n/a 

n/a 

n/a 

n/a 

Ho 

51080 

2 

6.78 

4.11 

0 

48.3 

10.0 

3 

31 

Ho 

51093 

2 

6.78 

3.84 

0 

49.2 

7.4 

2 

22 

Ho 

51095 

2 

6.78 

3.96 

0 

40.0 

10.0 

16 

5 

Ho 

60171 

0 

6.48 

3.73 

2.15 

n/a 

n/a 

n/a 

n/a 

Ho 

60558 

2 

6.78 

3.90 

0 

> 80 

8.3 

20 

23 

Ho 

60571 

2 

6.80 

4.28 

0 

> 80 

7.6 

8 

31 

Ho 
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Cluster ID 

Cat 

SED log(age) 

Log(mass) 

Av 

Max distance 

(pc) 

Min distance 

(pc) 

% max 

% min 

Best fit 

60582 

0 

6.60 

4.06 

3.48 

n/a 

n/a 

n/a 

n/a 

Ha 

60770 

1 

6.79 

4.09 

0 

73.9 

0 

10 

24 

Ho 

61053 

0 

6.70 

3.87 

3.08 

n/a 

n/a 

n/a 

n/a 

Ho 

61216 

2 

7.00 

4.13 

0.07 

> 80 

16.1 

21 

24 

Ho 

61310 

1 

6.92 

3.81 

0.32 

> 80 

4.5 

5 

11 

Ho 

61417 

2 

6.78 

4.57 

0 

> 80 

12.5 

10 

25 

Hq 


Table 1. Properties of the final clusters in the YMC sample of MSS, after removing poor sources. Cluster ID corresponds to the ID 
in the main Silva-Villa et al, 2014 catalogue, with the first digit indicating the field where the cluster is located. The ’cat’ indicates 
the category the cluster was sorted into: ’O’, ’1’ and ’2’ correspond to embedded, partially embedded and exposed respectively. Clusters 
indicated with a * are clusters that initially appear to be embedded from the H« morphology, though are partially embedded when you 
consider their SEDs, position on the colour-colour plot and their extinction. They seem to have cleared gas within our line of sight so 
they appear more exposed without visual inspection of the image. ’Max distance’ is the furthest radial distance cleared of gas by the 
cluster, ’Min distance’ is the minimum. ’% max’ and ’% min’ relate to the percentage of the circular region surrounding each cluster 
cleared to the maximum and minimum distances. Measurements of these properties were not obtained for embedded clusters. ’Best fit’ is 
the optimum SED fit for the cluster after checking each cluster’s individual SED, position on the colour-colour plot and Ho morphology. 
’Ho’ is a fit with U, B, V, Ho and I bands and ’No Ho’ includes U, B, V and I. 


at providing accurate ages for clusters younger than ~ 10 
Myr. Such inaccuracy means that any errors associated with 
studies of YMCs with ages in this range are likely to be 
significant. 

A novel method for providing a more accurate age, or 
for confirming previous age estimates has been developed by 


Gazak et al.|(2013), which uses the presence of red super¬ 


giants (RSGs) in the cluster as an absolute age indicator. 
By plotting the J-H colour for each of the clusters against 
log(age) (obtained via SED fitting), a clear distinction be¬ 
tween younger and older clusters exists in colour space. 
Younger clusters that are devoid of RSGs appear distinctly 
bluer, with values jumping to redder colours at older ages, 
immediately at the onset of RSGs. This is due to RSGs 
dominating the flux at red wavelengths, such as the near- 
IR. According to Gazak et al. (20131, any clusters with red 
J-H colours should be older than ~ 6 Myr, the age of the 
jump in the plot. One caveat for this method is that extinc¬ 
tion can also cause young clusters to have red J-H colours. 

For this study, the log(age) (based on our best fit, as 
quoted in Table was plotted against V-H colour (which 
was corrected for extinction) to find whether this jump is 
visible. The V band was chosen as no J band data was avail¬ 
able for WFC3 in Fields 3, 4, 5, 6 and 7. The plot is shown 
in Fig.[^ and despite the difference in wavelength, the same 
feature is produced with the jump in colour at the same age, 
~ 6 Myr. Almost all of the YMCs are where expected; older 
exposed clusters are redder in V-H, as they have RSGs and 
embedded younger clusters are the bluest. This indicated 
that our age estimates are fairly accurate, and at the values 
we would expect for these clusters. Though our SED fitted 
ages are rough estimates, this trend is still clearly observed, 
indicating the robustness of this technique. 


2 . 4.5 Wolf-Ray et features in clusters 

Wolf-Rayet (WR) stars are the later evolutionary stages of 
massive stars (typically 10-25 M©) and can be detected in 
clusters by broad emission lines in their spectra that are 
specific to these types of stars. The majority of WR stars 
are derived from 0-type stars and so have total lifetimes of 



Figure 8. V-H colour for all clusters in the M83 cluster sample, 
with the subsample of YMCs highlighted in colour with green 
stars, purple diamonds and orange circles corresponding to ex¬ 
posed clusters, partially exposed clusters and embedded clusters 
respectively. The jump from bluer, younger clusters to redder, 
older clusters due to the appearance of red supergiants is clear, 
with the exposed clusters at redder colours due to their slightly 
older age, and the youngest embedded clusters as the bluest. The 
model plotted over the data is calculated for the V-H band as per 
[Gazak et al.|p013| |. 


~ 5 Myr, 10 % of which is spent in a WR phase. Higher 
mass (25-30 M©) RSGs can also become WR stars. The 
characteristic emission lines arise from fast stellar winds and 
heavy mass loss from the star | Crowther|2007 1. 

The presence of such high mass, short-lived stars in clus¬ 
ters within our sample indicates that the clusters should be 
young. In older clusters, any WR stars would have already 
evolved past supernovae, given their ~ 5 Myr lifetimes, giv¬ 
ing us a refined age estimate of these clusters between 4-7 
Myr. [Bastian et al.j | |2006P previously used the same tech¬ 
nique to age date clusters in the Antennae galaxies. 


Hadfield et al. (2005 \ carried out a survey searching for 
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Properties of clusters with WR features 


Cluster ID 

Cat 

Log(Age) 

Log(Mass) 

10141 

1 

6.60 

4.03 

10438 

1 

6.85 

4.79 

10597 

2 

6.60 

4.72 

31015 

2 

6.60 

3.93 

50160 

1 

6.70 

3.70 

51080 

2 

6.78 

4.11 


Table 2. Clusters that coincide with identified WR regions in 
M83. ’Cat’ refers to classification of the cluster as embedded, 
partially embedded or exposed, as used throughout the paper. 
The ages and masses are the revised figures as described in § |2.4| 


WR stars in regions within M83, discovering 132 sources. We 
cross-referenced the coordinates of the regions they identi¬ 
fied as containing these objects with the coordinates of our 
cluster sample, examining images of the position of the fea¬ 
tures to conhrm that the features originated from the clus¬ 
ter. We found 5 regions that coincided with the position of 
clusters in the YMC sample. In addition to these 5 clusters, 
we also had optical spectroscopy of one further cluster in the 
sample that showed WR features. These clusters and their 
properties are shown in Table 

The low number of sources corresponding to clusters 
in our sample is mostly due to the lack of coverage of the 
original study. Their spectra did not uniformly cover the 
entire galaxy and so many clusters will have been missed. 
Additionally, these stars will be harder to detect in clusters 
due to being ’’washed out” by the continuum of the cluster, 
so many WR sources will be required for a detection. There 
are many young massive clusters that do not align with the 
WR sources, although it is not guaranteed that all clusters 
should contain enough WR stars to make detection possible. 
So the clusters without corresponding WR features may not 
be older, they just may not have sufficient WR stars to be 
detected. Therefore the WR feature method can only be 
used to confirm young ages for some clusters but not decide 
on the correct age. We believe the other methods discussed 
in this section should primarily be used for this purpose. 

Using this method indicates that 5 of the clusters with 
WR features have been correctly assigned younger ages. It 
also suggests that cluster 10438 should be younger. This 
cluster is already thought to have been given too old an age 
due to degeneracy in the model track, and this test agrees 
with that conclusion. The sample of 6 clusters used here 
show that WR features can be used to constrain ages of 
young clusters. 


3 RESULTS 

3.1 Statistical age results 

After comparison of the categorisations of the clusters with 
their age estimates, the results indicate that clusters are no 
longer embedded, and have removed their gas by ~ 4 Myr. 
The median ages for the embedded, partially embedded and 
exposed clusters were ~ 4 Myr, ~ 5 Myr and ~ 6 Myr, 
respectively. Therefore, the age at which clusters become 
exposed should be between 4 and 5 Myr. The maximum age 


of the clusters in the embedded sample is ~ 6 Myr, possi¬ 
bly indicating that the timescale for the clusters to progress 
from an embedded to exposed state is very short. Alterna¬ 
tively, this could just be an indication of the inaccuracies of 
SED htting at ages < 10 Myr, as discussed in § |2.4[ Older 
ages for the embedded and partially embedded clusters could 
arise from model degeneracies in SED fitting. These incor¬ 
rect ages may have skewed the median ages for embedded 
and partially embedded clusters, which should in fact be 
younger. This issue can be somewhat mitigated by deter¬ 
mining the correct fit, taking into account that embedded 
clusters should be younger than exposed clusters (therefore 
younger than ~ 6 Myr), as per the results of this paper. 

The majority of the clusters (53%) in the sample are ex¬ 
posed. If we include clusters that have removed any amount 
of gas from their surroundings i.e. partially embedded clus¬ 
ters as well, this means that 77% of the sample are no longer 
embedded. If we assume this sample is complete for M83 for 
clusters younger than 10 Myr (a 10 Myr age cut was ap¬ 
plied to the catalogue), and that the cluster formation rate 
is constant, this indicates that 7.7 Myr of the first 10 Myr 
of the life of young massive clusters are spent exposed to 
some extent. Alternatively, the first 2.3 Myr are spent in 
an embedded state, and clusters initially become exposed at 
this age, in agreement with the age indicated by the colour- 
colour plot discussed in § [2Z3 

A key assumption of the calculation, however is that all 
clusters are optically visible and enter our catalogue when 
they form. It can be argued that far younger clusters may be 
completely obscured and not be included, however a study 
by Whitmore & Zhang (20021 of clusters in the Antennae in¬ 
dicates that by comparing compact radio sources (extremely 
young clusters) with optical clusters, there is ~ 85% over¬ 
lap. This indicates that very few of these sources would be 
missed. Additionally, we searched the H band images, (as 
well as ground-based JHK band images in [Bastian et al.| 
(2014a l) and found no clear massive cluster progenitors that 
were not already included in our optical sample. This is also 
an indication that the embedded phase of the clusters is very 
short and that the clusters are very efficient at removing gas. 


3.2 Models of bubble expansion 


Dale et al. (20141 carried out simulations of cluster forma¬ 


tion and early evolution. Their simulations include the ef¬ 
fects of both ionisation and stellar winds from OB stars, 
although no radiation pressure has been included, which is 
likely to be relevant for very dense star-forming regions, e.g. 
[Fall et ^ ( |2010[ ). The simulations end 3 Myr after the for¬ 
mation of the first O-stars in each model cloud and some 
simulations already display bubbles formed by clusters with 
masses similar to the lower end of the range represented by 
our sample. This agrees with our results that these clusters 
will have removed their gas by 1-3 Myr. However, the bub¬ 
bles created in [Dale et al.| ( |2014[ ) are substantially smaller 
than those we observe. 

Having measured the maximum distance cleared by 
each of the exposed clusters, 77% of clusters had cleared 
some extent of their gas to radii greater than 80 pc. In con¬ 
trast, simulations of clusters with similar initial masses clear 
bubbles to a maximum possible distance of 40-50 pc in di¬ 
ameter. This indicates that the simulations do not reproduce 
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Figure 9. The plot of log mass of the cluster against the weighted 
average of the interior size of the bubble for exposed clusters. As 
shown, there appears to be no correlation between the two vari¬ 
ables, and may require a third variable to constrain a relationship. 


the extremely efficient process these clusters use to remove 
gas. It could mean that the process is related to radiation 
pressure, as this was not included in the simulations. Only 
a few of the simulations in Dale et al. (20141 are likely to 
be dense enough for radiation pressure to be important, but 
that may not be true of the clusters studied here. Alterna¬ 
tively, it could be an already known physical process that is 
far more effective than previously thought (such as stellar 
winds), or a completely new process. 

An important note is that these simulations end 3 Myr 
after the birth of the hrst O-stars. Should they run until 
6 Myr, the size of the bubbles in the simulations will in¬ 
crease and likely agree with our measurements. Therefore, 
we conclude that some of the simulations of gas removal 
are in reasonable agreement with what we find in the young 
clusters in M83. 


3.3 Cluster mass, distance cleared and cluster 
blowouts 

An interesting aspect of this study to take further, is whether 
or not there is a relationship between the mass of the cluster 
and the size of the bubble cleared. Using the measurements 
of the amount of space cleared by the cluster from Table 
we looked at the plot of mass of the cluster against the 
weighted average of the maximum and minimum distances 
cleared, as shown in Fig. There was no clear correlation 
between the amount cleared and the mass of the cluster, 
for clusters in our mass range of ~ 5000 — 6*10"^ M©. It is 
probable that if a relationship exists, it also requires a third 
variable to constrain it, such as position of the cluster in the 
galaxy (whether upstream or downstream of an arm, or in 
the inter-arm regions), or the cluster age. Including the clus¬ 
ter’s age in our plots still did not constrain any relationship, 
so we then decided to test the inclusion of galactic position. 

To test this theory, we visually identified the position of 
each cluster in the sample. This property was quantised into 
6 possible categories: in the bulge/inner bar region, in the 


outer bar region, at the end of the bar (which has increased 
levels of star formation), in the upstream regions of the arm 
just after gas enters the arm region and is shocked (lead¬ 
ing to star formation), in downstream regions at the arm’s 
outer edge, after the expected star formation, or outside of 
the arm overdensities. We plotted cluster mass against the 
weighted average distance of the cleared bubble as a func¬ 
tion of position in the galaxy and still no clear relationship 
could be determined. 

Including the mass of the gas outside of the exposed 
clusters could help constrain a possible relationship. How¬ 
ever, obtaining an accurate estimate of the amount of gas 
associated with the cluster would be very difficult using im¬ 
ages that rely on photometric bands highly susceptible to 
the continuum. Furthermore, the physical size of the shell 
would be very hard to determine, especially in more crowded 
environments. 

An additional aspect of the clusters that can be stud¬ 
ied is the presence of a ‘blowout’ in the gas morphology. 
We would define a blowout as free-streaming in one direc¬ 
tion within the gas cleared. Each of the exposed and par¬ 
tially embedded clusters were visually inspected to obtain an 
estimate of the number of clusters within the sample that 
contain a blowout region. We found that 18% of the en¬ 
tire cluster population had a noticeable blowout in the gas. 
As it would be impossible for embedded clusters to have a 
blowout (they have not yet strongly affected their environ¬ 
ment), when considering only partially embedded and ex¬ 
posed clusters, 24% of the population had a blowout. This 
indicates that this phenomenon should be fairly common 
for clusters of this mass. Pellegrini et al. (2012 I also studied 
blowouts in high mass (luminosity) HII regions in the LMC 
and found all regions were leaking photons. 


3.4 Implied constraints on age spreads within 
clusters 


Phenomena observed in intermediate age (1-2 Gyr) clusters 
in the LMC and SMC, such as double, or extended main se¬ 
quence turn offs, e.g. [Mackey fc Broby Nielsen|(|2007[),|Goud-| 


frooij et al. (20091, Girardi et al. {20131, and in globular clus¬ 


ters such as split main sequences, e.g. Piotto et al. (20051, 
[Milone et al.| ( |2013| ) , and chemical spreads have widely been 
explained using age spreads, or multiple bursts of star for¬ 
mation. The extent of the age spreads can be up to a few 
hundred Myr (e.g. [Milone et aL| (|2009P), thoug h this has 
been strongly contested by Bastian & Silva-Villa (20131 and 


Niederhofer et al. (20151, who have found no evidence of 


age spreads in YMCs with resolved photometry. Addition¬ 
ally, stellar rotation has been shown to potentially exhibit 


the same effects seen in the main sequence turn offs {Bastian 
[fc de Mink|2009[ | . [Cabrera-Ziri et ^ ( [2014[ ) also showed that 
the spectra of NGC 34 can only be reliably ht with a single 
stellar population, contesting the theory of multiple bursts 
of star formation. 

The debate of the origin of multiple stellar populations 
continues, and our results can be used to apply some con¬ 
straints to current claims of age spreads. The clusters stud¬ 
ied here are far younger than those with the proposed age 
spreads, and they have already become exposed and free of 
gas at 4 Myr. The lack of ionised gas within the cluster at 
this time means that there can be no star formation until 
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more gas has been accreted, although the mechanism used 
to accrete this gas is still unclear. 

As shown in Bastian et al. (2014bI, this is also true 
of very young clusters, of approximately the same age but 
larger masses than the clusters studied in this paper. In that 
study it was shown that the gas was still moving away from 
the clusters at a speed of ~ 40 km/s, and therefore the pris¬ 
tine gas required to produce the second generation of stars 
would need to be accreted from other sources, unrealistic for 
the amount required. Our clusters show very similar results. 
The removed gas should be receding from the cluster at a 
speed of ~ 20 km/s, and clusters of this mass would not be 
able to efficiently re-accrete this gas in the proposed time 
between the two generations forming. It is highly unlikely 
the cluster would be able to accrete this required material 
from elsewhere in the galaxy on this timescale (e.g. |Bastian| 


et al. (2014a l) also, especially considering the amount of gas 


required would be on par with the existing cluster mass for a 
reasonable star formation efficiency. This indicates that the 
current required age spreads within clusters are unrealistic 
for young massive clusters in this sample. 


4 DISCUSSION 


The results of this paper have important consequences for 
one of the current globular cluster formation scenarios, the 


Fast Rotating Massive Star Scenario (FRMS) (Krause et al. 
2013). In this scenario, the hrst generation of stars form. 


which then become mass segregated leaving the high mass 
stars in the centre. Bubbles of gas surround each high mass 
star in narrow filaments. Decretion discs form around the 
high mass stars and ejected material mixes with the pristine 
gas from the filaments to form the second generation of stars 
within the disks. When the high mass stars have evolved to 
the supernova stage, the leftover pristine gas from the fila¬ 
mentary structures is arranged in a larger bubble surround¬ 
ing all stars in the cluster, so that the cluster appears to be 
embedded in gas. 

The young clusters must remain embedded in their na¬ 
tal gas for 25-30 Myr in the FRMS scenario, however the 
clusters in this study have shown that this is not the case, 
and they have in fact removed their excess gas by < 4 Myr. 
This clearly indicates that the FRMS scenario, in its current 
form, is not feasible and clusters are far more efficient at re¬ 
moving gas than expected. One clear difference is the mass 
of the cluster that has formed. The FRMS scenario discusses 
clusters of mass > 10® Mq, whereas the most massiv e cluster 
in this work is ~ 10"^ M©. However, the results of Bastian 


et al. (2014b I suggest that the effect of clusters clearing out 


gas is scalable between cluster masses. They present similar 
results for clusters of masses ~ 10® M©, on par with those 
described in the FRMS scenario. The same effect is seen 
in these clusters; gas is cleared at the same age as for the 
lower mass clusters. One difference is the size of the bubble 
cleared. Higher mass clusters appear to create larger bub¬ 
bles on the same timescales, probably due to larger feedback 
acting to remove the gas. 

The mechanism acting within the clusters to remove the 
gas is much more significant than previously assumed, shown 
by the clusters ability to clear gas so early in their evolu¬ 
tion. Supernovae have been considered as one of the main 


contributors to the removal of gas. It is feasible to include su¬ 
pernovae for systems at approximately solar metallicity, like 
M83, however this only holds if the cluster contains very high 
mass stars. For the first supernovae to occur at 2-3 Myr after 
cluster formation (the age at which the clusters seem to start 
removing gas) the IMF must be well sampled. A reasonable 
estimate for the most massive star in clusters of ~ 5000 M© 
is ~ 40 M©, which would not undergo a supernova phase 
until ~ 8 Myr, and so will not have made an impact on the 
removal of gas. Therefore for supernovae to be considered, 
stars much greater than this mass (likely at least one star 
of ~ 120 Mq) must exist in each of the clusters, which is 
unlikely to be the case for all clusters. At sub-solar metallic- 
ities, these very high mass stars would collapse straight into 
a black hole phase, but supernovae could at least make a 
contribution to the gas removal at solar metallicity, though 
they are very unlikely to be the sole mechanism |Heger et al.| 
( 20 ^. 


5 CONCLUSIONS 

Here we summarise our main conclusions: 

• The young massive clusters in our sample are all ex¬ 
posed and free of gas at ~ 4 Myr. The process of removing 
gas seems to begin at 2-3 Myr, as indicated by the position 
of exposed clusters on in colour-colour space and statistical 
considerations of the relative population sizes of exposed or 
partially embedded clusters as opposed to embedded clus¬ 
ters. 

• The mechanism for the removal of gas from clusters 
is far more efficient than previously thought. Supernovae 
could play a role but cannot be the sole contributors. Ra¬ 
diation pressure could contribute, as well as strong stellar 
winds from stars at advanced evolutionary stages, such as 
Wolf-Rayet stars, which have been identihed in several of 
the YMCs in the sample. 

• Current models of gas clearance from clusters could 
potentially produce similarly sized bubbles when compared 
with our images, though radiation pressure was not included 
and could be a key contributor to the removal of gas. 

• The lack of gas within these clusters provides con¬ 
straints for the age spreads that have been proposed to ex¬ 
ist. A spread of a few hundred Myr is required, however 
these clusters have removed their gas by ~ 4 Myr, and so 
cannot form a new generation of stars, unless this gas is re¬ 
accreted, or new pristine gas is obtained. The amount of gas 
that would need to be accreted over this short time frame 
is unlikely, especially considering the cluster’s original gas is 
pushed away to large distances, out of reach for re-accretion. 
It is also unlikely that gas from other sources could be ac¬ 
creted in the space of a few hundred Myr. This also has con¬ 
sequences for the FRMS scenario which currently requires 
clusters to remain embedded in natal gas for 25-30 Myr in 
order to form a second generation of stars is infeasible, as 
these clusters have removed gas by ~ 4 Myr. 

• Age fitting using SEDs is often unreliable for young 
clusters due to degeneracies in the models and aperture ef¬ 
fects. Constraints on age can be made by using the near- 
infrared bands, such as H (F160w), where the characteristic 
jump at ~ 5 — 6 Myr from blue to red colours due to RSGs 
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can provide an upper limit. The presence of Wolf-Rayet fea¬ 
tures also puts a constraint on age of the cluster of 4-7 Myr. 
Using Ha in the fit for the age of a young cluster and as an 
indicator of age is very unreliable as it is not directly linked 
to the stars, but rather line emission and the continuum can 
have an effect. A better indicator could be the UV. 

• No relationship has been found between the mass of the 
cluster and size of the bubble cleared for clusters of mass 
between 5000 and 6a:10^ M©. If a relationship does exist, it 
likely requires more variables to constrain the correlation. 

• Approximately 25 % of the clusters were identified as 
having ’blowouts’. This indicates that it is a fairly common 
phenomenon in clusters. 
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